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ABSTRACT 

There is a well known correlation between the mass and metallicity of star-forming galaxies. Because 
mass is correlated with luminosity, this relation is often exploited, when spectroscopy is not available, 
to estimate galaxy mctallicities based on single band photometry. However, we show that galaxy color 
is typically more effective than luminosity as a predictor of metallicity. This is a consequence of the 
correlation between color and the galaxy mass-to-light ratio and the recently discovered correlation 
between star formation rate (SFR) and residuals from the mass-metallicity relation. Using Sloan 
Digital Sky Survey spectroscopy of 148,021 nearby galaxies, wc derive "LZC relations," empirical 
relations between metallicity (in nine common strong line diagnostics), luminosity, and color (in three 
filter pairs). We show that these relations allow photometric metallicity estimates, based on luminosity 
and a single optical color, that are ~ 40% more precise than those made based on luminosity alone; 
galaxy metallicity can be estimated to within ^ 0.06 — 0.1 dex of the spectroscopically-derived value 
depending on the metallicity diagnostic used. Including color information in metallicity estimates 
also reduces systematic biases for populations skewed toward high or low SFR environments, as we 
illustrate using the host galaxy of the supernova SN 2010ay. This new tool will lend more statistical 
power to studies of galaxy populations, such as supernova and gamma-ray burst host environments, 
in ongoing and future wide field imaging surveys. 

Subject headings: galaxies: abundances — galaxies: photometry — galaxies: ISM 



1. INTRODUCTION 

The gas-phase metallicity of galaxies, as mea- 
sured from their nebular emissi on spectrum, is cor- 
related with galaxy lu minosity (jLequeux et al.l 119791 : 
iGarnett fc Shields! Il987| ). This relation has been used 
as a key observational tool in the study of populations 
such a s supernova host ga laxies fe.g. lPrantzos fc Boissieii 
l2003t lArcavi et al.l I2O1O0 . where gas-phase metallicity 
is an important proxy for the properties of their short 
lived progenitor stars. However, using Sloan Digital 
Sky Survey (S PSS) imaging a r id sp ectroscopy for ^ 
53,000 galaxies. iTremonti et all (|2004[ ) showed that the 
luminosity-metallicity relation has a large intrinsic scat- 
ter of cr = 0.16 dex (50%), in terms of metallicity residu- 
als, which limits the utility of this relation as an effective 
indicator of metallicity. 

There are two primary causes for the scatter in the 
luminosity-metallicity relation. First, while the scatter in 
the mass-m ctallicity relation is fairly small (cr = 0.10 dex, 
ITremonti et al. .2004ii) . luminosity is not a perfect proxy 
for mass. The mass-to- light ratio of galaxies is highly cor- 
related with galaxy color, such tha t redder galaxies at a 
fixed luminosity are mor e massive (|Bell fc de ,Tongil200ll : 
iKauffmann et al.ll2003bl ). Second, a more fundamental 
relation has been uncovered between mass (M), metallic- 
ity (Z), and star forma.tion r ate (SFR) ()Lara-L6pez et ahl 
120101: lIMannucci et al.|[2010D . This "fundamental plane" 
or "Fundamental Metallicity Relation" has remarkably 
small residual scatter (cr = 0.05 dex), indicating that 
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variations in SFR are responsible for much of the scatter 
in the mass-metallicity relation. The existence of this 
fundamental plane is a valuable constraint for models of 
galaxy evolution, and likely an expressio n of galaxy out- 
flows, infall, and/or dow nsizing (Mannucci et aLl boiOl: 
IPeeples fc Shankadl20ll . To improve the precision of 
photometric metallicity estimates, a readily accessible 
observable must be used to break the degeneracy between 
luminosity, mass, and SFR. 

In this paper, we show that the addition of color infor- 
mation significantly decreases the scatter in photometric 
metallicity estimates. We derive the optimal projection 
of the fundamental plane for star-forming galaxies, in 
terms of the observable properties luminosity and color, 
that we call the LZC relation. Wc describe the sam- 
ple of SDSS galaxies we use to study these correlations 
and methods for spectroscopic metallicity estimation in 
Section [S] In Section [3l we derive analytic expressions 
for the LZC as expressed in a variety of different filter 
sets and metallicity diagnostics. Finally, in SectionjH we 
describe how specific observational studies may benefit 
from the LZC relation in making precise metallicity de- 
terminations from imaging available from wide field sky 
surveys. 

2. GALAXY SAMPLE 

We used the MPA/JHU catalog of 927,552 star- 
formi ng galaxies from the SDSS-DR7 (lAbazajian et al.l 
l2009f ). The catalog includes synthesized rest- frame pho- 
tometry corrected for Galactic extinction, emission line 
fluxes, stellar masses (based on SED fitting to ugriz pho- 
tometry) , and star formation rates fo r each galaxy, as de- 
scribed in IKauffmann et all ()2003bl ) ; iBrinchmann et aD 

The MPA/JHU catalog is available at 
|http://www. mpa-garching.mpg.de/SDSS] 
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TABLE 1 



Metallicity Diagnostics Used 


Name 


Line Ratios 


Source 


M91 


R23 


McGaueh 1991'' 


POl 


P 


Pilvuein 2001 


KD02 


N202 


Kewlev & DoDita 2002'' 


KK04 


R23 


Kobulnickv & Kcwlc^^ 2004 


PP04 


N2,03N2 


Jcttini & Pagcl 2004 


T04 


model'^ 


I'remonti ct al. 2004 


PT05 


P 


Pilvuein & tliuan 2005 


PVT 


P 


Pilvuein et al. 2010'' 



We have adopted the rev ised prescription suggested 
bv lKobulnickv et al.l 119991) . 

FoUowing lKewlev &: Ellisoiil II2008I '). we adopt the av- 
erage of the M91 and KK04 for the lower branch solu- 
tion. 

The T04 metallicities are estimated based on simulta- 
neous fits of all major emission lines to photoionization 
models and are provided in the MPA/JHU catalog 

We use the "ONS" solution, which includes a de- 
pendence on the S II flux, for the conditions log(N2) > 
-0.1 and log(N2/S2) > -0.25 (which is true for ~ 98% 
of the SDSS galaxies). 

(|2004[) : ISalim et all (|2007[ ). We correct line fluxes for 
line-of-sight extinction (including both Galactic and 
intrinsic reddening) from the Balmer flux decrement. 
We assume Fhq/Fh/3 = 2.85 (corresponding to T = 
10, OOP K and Up = IQ^ cm~ ^ for Case B recombination; 
Osterbrock & Fcrland '2006") and the extinction curve of 
Cardelh et aTTiiigsa '). assuming Ry = 3.1. 

We perform preliminary s ample cuts on the cat alog 
following the prescription of iMannucci et all ()2010[) . as 
follows. We require 0.07 < z < 0.3, (S/N)h„ > 25, 
Fy[cy/Fy[fj > 2.5, Ay < 2 . 5 mag, and we reject AGN 
following iKauffmann et all (|2003a[ ). In total, following 
these cuts, we consider 148,021 galaxies. 

For each galaxy, we compute oxygen abundance as a 
proxy for metallicity using a variety of strong line di- 
ag nostics that are w idely used in the literature (Table (TJ 
see lLopez-Sanchez et al. 20i3for a recent review). First, 
we employ several diagnostics relying on the R23 ra- 
tio, which depends on the fluxes of O II AA3726, 3729, 
O III A4959 and A5007, and H^. The R23 diagnos- 
tic su ffers from a degeneracy (see e.g. iKewlev fc Dopital 
I2002D that we break using either the N202 (N II A6584 
to O II A3727) or N2 (N II A6584 to Ha) ratios, fol- 
lowing the authors' prescriptions. Next we employ di- 
agnostics depending on N2, N202, and 03N2, the flux 
ratio of O III A5007 and N II A6584. Finally, we em- 
ploy diagnostics based on the ionization parameter, P, 
the ratio of O III A4959 and A5007 to R23. It is nec- 
essary to use multiple diagnostics because they exhibit 
well known systematic discrepancies, which arc particu- 
larly strong between those diagnostics calibrated empiri- 
cally and those calibrate d against photoionization models 
(IKewlev fc EllisonI 120081) . 

Figure [1] demonstrates two correlations in the galaxy 
sample. First, it shows the well known luminosity- 
metallicity correlation, where more luminous galaxies 
(shown in Mr) typically have higher metallicities. How- 
ever, there is significant scatter in this relation, with a 
standard deviation of <to = 0.14 dex in the residuals 
in metallicity for r band, as shown. Second, the fig- 
ure demonstrates that there is a correlation between the 
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Fig. 1. — A luminosity-metallicity (r-band, T04) plot showing a 
random subset of 10, 000 SDSS galaxies, color coded hj g — r color. 
The solid lines show the median luminosity-metallicity relation for 
the galaxies divided into 10 equal-sized bins in g — r color. 

residual in metallicity (the offset from the luminosity- 
metallicity relation) and galaxy color. 

3. LUMINOSITY-METALLICITY-COLOR RELATION 

Following IMannucci et al.l ()2010[) . we project the LZC 
relation onto an axis /i with components of color and 
luminosity: 

fi = Mi - a X [rrij - rrik) (1) 

12 + log(0/H) = po + Pi M + P2 + P3 (2) 

where i , j, A; are choices of filters and pi are parameters of 
a third order polynomial. For each combination of metal- 
licity diagnostic, luminosity band, and color, we deter- 
mine the optimal projection of the galaxy LZC relation 
by sampling from a grid of a parameters and selecting 
the value that minimizes the variance in the residuals of 
metallicity. We report the optimal value of a and corre- 
sponding best fit LZC parameters p in Table [H 

We find a correlation between the luminosity- 
metallicity relation and color which varies in strength 
depending on the choice of metallicity diagnostic and 
filter set. Figure [2] demonstrates this optimization for 
one filter set and three choices of metallicity diagnos- 
tic: T04, KD02, and PP04 03N2. For PP04 03N2, the 
scatter in the metallicity residuals of the LZC relation 
is az = 0.10 dex (^ 25%), as compared to the LZ re- 
lation (a — 0), az.o = 0.13 dex ('^ 35%). The decrease 
in residual scatter is similar in other diagnostics, ranging 
from az/<7z,o = 0.68 — 0.87, with the exceptions of POl 
and PT05 (discussed below). Note that, regardless of the 
choice of diagnostic or filters, the residual scatter is lower 
for asymptotically high values of a than it is for a = 0. 
This implies that, in general, color is more effective than 
luminosity as a predictor of metallicity. 

In contrast, the residual scatter achieved by 
IMannucci et al.l ()2010f ) in terms of the optimal projec- 
tion of the physical parameters stellar mass (Af*) and 
SFR, /Lio.32 = log(M*)-0.321og(SFR), was only 0.05 dex. 
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Fig. 2. — The optimal projection of the LZC relation for Mr, g — r color, and three different metallicity diagnostics (T04,KD02, and 
PP04 03N2). The projected LZC relation is shown for the optimal value of a, with the best fit third order polynomial LZC relation 
in black, extending ov er the calibrated range (2 — 98th percentile). The color coding shows the optical physical parameter {(132) from 
IMannucci et al.1 I I2010I) . The insets display the dispersion in the residuals in metallicity from the LZC as a function of the color- weighting 
parameter a, with the optimal value marked by the dashed line. 



However, estimation of /io.32 is based on full ugriz imag- 
ing and R ^ 2000 optical spectroscopy, while the LZC 
relation relies on imaging in just two bands and a red- 
shift estimate. The coloring in Figure [5] illustrates that 
the optimal projection of the photometric properties is 
highly correlated with /io.32, with Pearson correlation co- 
efficient p ^ 0.8 for all metallicity diagnostics and filters. 

For two of the metallicity diagnostics we consider (POl 
and PT05, both based on P), the strength of the LZC is 
negligible (ctz/ctz.o ~ 1). This could imply that the in- 
formation expressed by galaxy color is already included 
in these metallicity diagnostics, such that they produce 
more physically realistic metallicity estimates than the 
other diagnostics. However, we find the opposite sce- 
nario to be more likely. These metallicity diagnostics 
have well known discrepancies with other strong line di- 
agnos tics, yielding the M-Z re lation essentially unobserv- 
able (iKewiev fc EUisonl 120081 ). Similarly, we find that 
they produce extreme scatter when used to calibrate the 
luminosity-metallicity relation {<tz.o = 0.19 dex for POl, 
Mr, {r — i)), which is not improved by the incorporation 
of color information. However, the LZC relation cali- 
brated with the PVT diagnostic (also based on P) does 
exhibit significant decrease in scatter {(Jz/<Jz,q ~ 0.85 
depending on filter set). 

As expected, choosing a bluer filter set produces the 
greatest improvement in uz because it better reflects the 
size of the young stellar population. Taking the median 
across all metallicity diagnostics and using Mi, we find 
c^z/crz,o = [0.81,0.82,0.88] for [g ~ r,g - i,r - i]. The 
choice of filter for L determination, among gri, has a less 
significant effect. Using the g — i color, we find uz /<^z = 
[0.78,0.80,0.82] for [5, r,i] 



4. APPLICATIONS 

The calibrated LZC relation we present could ben- 
efit several disciplines, allowing for precise and accu- 
rate metallicity estimates for galaxies based on photom- 
etry alone. Spectroscopic metallicity measurement de- 
mands considerably mor e observational reso urces, while 
full SED modeling (e.g. iPacifici et all [20121 ) IS compu- 
tationally demanding. In contrast, the LZC relation 
can be employed using existing multi-band photome- 
try from wide-field surveys such as SDSS or newly ac- 
quired, targeted observations, so long as an estimate for 
the redshift of each galaxy is available. Using ugriz 
photometry, photometric redshift estimates for galax- 
ies ca n be achieved with scatter 5z ~ 0.02 ()Ball et al.l 
120081). while i-iT-correctio ns can be determined to < 20% 
( Blanton fc RoweisI [20071 ). 

To test the uncertainty in metallicity introduced by 
use of photometric redshifts (photo-z), we recompute 
metallicities for the subset of 70, 000 galaxies in the 
MPA-JHU catalog with photo-z estimates in the SDSS- 
DR8. We use the fcd-tre e nearest neighbor fit photo-z 
estimates, as described in lCsabai et al.l ()2007D . We com- 
pute j^-correctioii s using the analytic prescriptions of 
iChilingarian et all ()2010l ). using both the spectroscopic 
add photometric redshifts for each galaxy. We apply the 
LZC relation as calibrated for the PP04 03N2 metal- 
licity diagnostic using the Mg luminosities and g — r 
colors. The resulting distribution of metallicity resid- 
uals {5Z) for the spectroscopic versus photometric red- 
shifts suggests that there is no systematic bias introduced 
by photo-z (median 5Z = 0.000 dex). The typical un- 
certainty added to the metallicity estimates is negligi- 
ble (standard deviation 5Z = 0.008 dex) and therefore 
the photometric metallicity estimate is dominated by the 
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TABLE 2 

Parameters of the LZC relations 



L Color a P3 x 10^ P2 X 10'' pi x 10'' po o^z,o C"z o"z/o"z,o A* range'' 



PP04 03N2 



9 


9 - 


r 


28.7 


0.80 


-1.02 


-122.93 


6.01 


0.130 


0.088 


0.678 


-43.3, 


-26.5 


9 


9 - 


i 


11.8 


10.00 


5.13 


16.54 


7.33 


0.130 


0.088 


0.677 


-34.0, 


-22.4 


9 


r — 


i 


18.1 


2.88 


0.81 


-53.50 


7.37 


0.130 


0.099 


0.761 


-26.9, 


-19.1 


r 


9 - 


r 


27.7 


0.80 


-1.02 


-122.93 


6.01 


0.127 


0.088 


0.696 


-43.3, 


-26.5 


r 


9 - 


i 


10.9 


11.97 


6.42 


42.01 


7.47 


0.127 


0.088 


0.691 


-33.7, 


-22.4 


r 


r — 


i 


15.9 


4.15 


1.10 


-63.26 


7.14 


0.127 


0.098 


0.773 


-26.8, 


-19.4 


i 


9 ~ 


r 


24.2 


2.47 


0.14 


-97.70 


6.21 


0.123 


0.088 


0.712 


-40.8, 


-25.6 


i 


9 ~ 


i 


10.8 


10.00 


5.13 


16.56 


7.33 


0.123 


0.088 


0.714 


-34.0, 


-22.4 


i 


r — 


i 


14.9 


4.15 


1.10 


-63.26 


7.14 


0.123 


0.098 


0.795 


-26.8, 


-19.4 



Note. — Parameters of the LZC relation defined in Equation [T] L is the photometric band 
of the luminosity, a is the optimal value of the color weighting factor to minimize the scat- 
ter in metallicity, and pi are the parameters of the best-fit LZC polynomial for the correspond- 
ing value of a. The value az/crzfi expresses the reduction in the scatter relative to the LZ re- 
lation (without color term). Tabic [2] is published in its entirety in the electronic edition and 
at https://www.cfa.harvard.edu/-nsaiiders/papers/LZC/summary.htm, including parameters for all 
metallicity diagnostics. A portion is shown here tor guidance regarding its form and content. 
^ The range of the optimal projection of luminosity and color, over which the diagnostic is calibrated. 
This range is defined by the 2nd and 98th percentiles in n of the SDSS galaxies in the calibration sample. 



s catter in the LZC r elatio n. 

iLara-Lopez et al.l (|2010D concluded that there is no 
detectable evolution in the fundamental pl ane for star- 
forming galaxies out to z ^ 3.5. Similarly, iCresci et al.l 
(|20T1) inspect - 300 zCOSMOS galaxies and find no 
evolution in the Fundamental Metallicity Relation to 
z ^ 0.8. This would imply that observers who wish 
to apply the LZC relation to imaging of galaxies in 
this reds hift regime need only accoun t for i^T-corrections. 
However. iPerez-Montero et al.l (j2012| ) investigate a larger 
sample of ~ 5000 zCOSMOS galaxies to z ~ 1.3 and re- 
port evolution of the SFR-corrected mass-metallicity re- 
lation starting a,t z > 0.4. This suggests that additional 
studies are needed to calibrate the LZC for redshifts 
z > 0.4. 

Studies of supernova (SN) host galaxies c an support in- 
feren c es into pr ogenitor st ar populations ("Modia z et alJ 
[2008t iPrieto et al. 2008t iSanders et all 2012b, e.g.), 
with some studies relying on photometry rather than 
spectroscopy to mea sure t he host galaxy met a llicity 
(iPrantzos fc Boissieri l2003t iBoissier fc Prantzqd 1200^ 
lArcavi et al.ll2010l e.g.). However. ISanders et al.l ()2012bl ) 
have shown that the statistical uncertainty associ- 
ated with metallicity estimates based on the galaxy 
luminosity-metallicity relation is a significant barrier to 
detecting subtle differences in metallicity among SN pop- 
ulations. Moreover, because some SNe strongly prefer 
blue galaxies (high SFR envi ronments; iLevesque et al.l 
120101 : iKellv fc Kirshneil l2011| ). photometric metalHcity 
estimates will be biased if color information is not in- 
corporated. Using the LZC relation will effectively re- 
move this bias, and significantly reduce the uncertainty 
in metallicity measurements. The additional uncertainty 
introduced by photo-z should be minor, as SN host galax- 
ies studies are almost exclu sively done in the z < 0.15 
regime (see compilation in ISanders et al.l l2012b[ ). and 
spectroscopic redshift estimates are often available from 
the SN spectroscopy. We note that the LZC relation pre- 
dicts the "global" or nuclear galaxy metallicity, as probed 
by the 3" SDSS spectroscopic fibers, and significant off- 
sets may exist from the SN host environment metallic- 



ity due to metallicity gradients in galaxies. However, 
these metallicity offs ets are typically small (< 0.1 dex; 
ISanders et al.ll20l"2B l. and the intrinsic scatter in galaxy 
metallicity profiles limits the abilit y to isolate the explo- 
sion site even with spectroscopy (|Rosolowskv fc SimonI 
HOOS; Sanders et al. 2012'a[). 

Similarly, galaxy metallicity measurements are key 
to the discussion of the progenitor properties of 
long-duration ga,mma-ray bursts (L GRBS, see e.g. 
iFvnbo et al.l[200l[SMld"ers et al.ll2012d ). While the host 
environments of LG RBs typically fall bel ow the mass- 
metallicity relation (jLevesQue et al.|[2010[ ). it has been 
shown that these h ost galaxies do follow t he fundamen- 
tal plane relation (iMannucci et afl 120111 but see also 
IKocevski fc Westll201in ." Because LGRBs are frequently 
discovered at high redshift (z > 1), the LZC relation 
could be used to derive metallicity estimates at consid- 
erably lower expense than deep NIR spectroscopy (e.g. 
iMaiolino et al.ll2008D . However, in order to study the 
extremely lowest metallicity environments preferred by 
LGRBs {Z < 0.3 Zq), additional calibration is needed 
to extend the LZC relation beyond the range probed by 

xy sarn i: 

IMannucci et alfeoilt) . 

As a usage example, we apply the LZC relation to the 
unusua l host galaxy of the SN 2010ay. In ISanders et alJ 
(|2012cl ). we report that this host galaxy is a 2a out- 
lier from the luminosity-metallicity relation. The median 
and 1(T metallicity interval for SDSS galaxies with lumi- 
nosity similar to this host galaxy {Mb = —18.3 mag) 
is 12 + log(0/H) = 8.93 ± 0.17 (T04). This is a fac- 
tor of > 2 greater than the spectroscopically-measured 
T04 metallicity of 12 + log(0/H) = 8.58. The discrep- 
ancy is due to the extremely low mass-to-light ratio and 
high SFR of the host galaxy. Using the LZC relation 
for Mr and r — i color, we find a T04 metallicity of 
12 + log(0/H) = 8.72 ± 0.09. This agrees with the spec- 
troscopically measured value to within ^ 40% and has a 
significantly lower associated uncertainty. However, the 
LZC relation cannot be applied in all filter combinations 
because the host lies outside the calibrated range for fi, 



the SDSS galaxy sam ple, 0.4 Zq < Z < 1.3 Zq (see e.g. 
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illustrating the importance of extending the calibration 
presented here to lower-metallicity host galaxies not well- 
represented in the SDSS spectroscopic sample. 

Finally, we suggest that the next generation of wide 
field, multi-band, photometric surveys could use the 
LZC relation to characterize the metallicity distribu- 
tion of galaxies in the local universe, and its evolution 
with redshift. The Pan-STARRSl (PSl) survey is al- 
ready operating, and will provide gpi rpi ipi zpi ypi 
photometry for ~ 2 x 10® galaxies over 3/4 of the sky 
(jSagha et al.l l2012| ) . Simultaneousl y, the Baryon Oscil- 
lation Spectroscopic Survey (BOSS iDawson et al.ll2012l ) 
is providing spectroscopic redshifts to ~ 1.5 x 10^ galax- 
ies to z ^ 0.7, but with emission line S/N typically too 
low for spectroscopic metallicity analysis. In the future, 
the Large Synoptic Survey Telescope (LSST) will pro- 
vide uqrizy photometry of ^ 10^*^ gal axies to z 6 
(|LSST Science Collaboration et al.l 120091 ). With the ad- 
vent of such datasets, the LZC relation may play an 



important role in defining the metallicity distribution of 
galaxies that has emerged from the cosmic evolution of 
star formation and galaxy mass. 
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